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DUKHIN AND PECLET NUMBERS IN THE 
PRESENCE OF SLIP 



ELECTRO-OSMOTIC VELOCITY IN 
EIGENDIRECTIONS 
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Surface conduction in presence of slip is characterized 
by full Dukhin number, which is given by 
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where m = 2e{kBT / ze)'^ / {i]D), and D is the ion diffusiv- 
ity, which is assumed equal for both types of ions. The 
first term in ([T]) is the Dukhin number for the no-slip 
surface, Duh=o, while the second one, I?uti, is due to hy- 
drodynamic slip. In the Debye-Hiickel model 



< 1. 
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This restriction ([2]) allows the simplification: 
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The parameter m k, 100/z^, and kL ^ 1 since EDL 
is thin. Whence Dub=o ~ \/{z^kL) <C 1 provided the 
potential is low. For a "slippery part" of Du we evaluate 
Dub « 0.56/^ if eC/{kBT) « 0.1, and Dub « 5 • lO'^b/L 
for e(/{kBT) « 0.01. Therefore for increasing surface 
charge (potential) and b/L, the conductivity of the diffuse 
layer can become comparable to the bulk, and surface 
condition must be considered. 

The Peclet number, Pe = UL/D , in presence of slip 
can be evaluated as 
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Typically, electroosmotic velocity is of order is of order 
micrometers per second for no-slip surfaces ■ For nano- 
scale patterns L < 1 ^m and typical ion diffusivities D « 
10~^ cm^/s this gives Peb=o < 0.01 <C 1. The slip implies 
a correction factor (1 -|- 6k) , which suggests that the 
convective ion transport can safely be neglected only for 
bK < 10. Larger values of 6k should relax this standard 
approximation of small Pe. 



Longitudinal stripes.- In this configuration only 
X— velocity component remains, and the Stokes equation 
takes the form 



(Sy +dl)u = EK^i^Et 



(6) 



We expand surface charge density in a Fourier series, and 
the potential is then 

V(2/,z) = Me-« + V^e-«"^cosA„z, (7) 
en sin 



where ^„ = k"^ + A^ , A„ = 2n7r/i , {q) = qi(f)i + q2(j)2 
is the mean surface charge, and 
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Qn = sm — — . (8) 
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The general solution to ([6]) for u(y, z) has the form 

°° eE 
u{y, z) = C/ll + y Une-^-y cos A„z + — (9) 

71=1 ' 

where C/„ and fj" are determined by the slip boundary 
conditions. Imposing them on ([S]) in a thin EDL limit 
yields a dual series 
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which can be solved exactly by using a technique [3| to 
obtain the thin-EDL electro-osmotic velocity: 
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Transverse stripes. -Although an external pressure 
gradient is equal to zero, local pressure variations con- 
tribute into a non-zero term Vp in the Stokes equation, 
so that the flow is essentially two-dimensional. We first 
introduce a stream function f{x,y) 



dyf = It, dxf 
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which obeys inhomogeneous biharmonic equation: 



These dual series can be solved exactly to obtain 
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Here u and v are x and y velocity components, corre- 
spondingly. The general periodic solution to has the 
form 
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Here the potential ^{x, y) has exactly the form ([7]) with 
z replaced by x. The dual series problem in a i/iin iS.DL 
limit can be written as 
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C/^ + a„ cos A„a 



5/2 < \x\ < L/2, (17) 



where 



We emphasize that a comparison of Eqs. ([T2)) and ([18]) 
indicates that the EO flow is generally anisotropic, so 
that our results do not support an earlier conclusion [3| 
that the electro-osmotic mobility tensor is isotropic in 
the thin EDL limit. This inconsistensy Q (due to an 
erroneous expression for a transverse electro-osmotic ve- 
locity, where factor of 2 was lost) has been corrected for 
a case 62 = oo in Q . 
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We give a general theoretical description of electro-osmotic flow at striped super-hydrophobic 
surfaces in a thin double layer limit, and derive a relation between the electro-osmotic mobility 
and hydrodynamic slip-length tensors. Our analysis demonstrates that electro-osmotic flow shows 
a very rich behavior controlled by slip length and charge at the gas sectors. In case of uncharged 
liquid-gas interface, the flow is the same or inhibited relative to flow in homogeneous channel with 
zero interfacial slip. By contrast, it can be amplifled by several orders of magnitude provided slip 
regions are uniformly charged. When gas and solid regions are oppositely charged, we predict a flow 
reversal, which suggests a possibility of huge electro-osmotic slip even for electro-neutral surfaces. 
On the basis of these observations we suggest strategies for practical microfluidic mixing devices. 
These results provide a framework for the rational design of super-hydrophobic surfaces. 
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Introduction.- Electro-osmotic (EO) "plug" flows 
are established when an electric field forces the diffuse 
ionic cloud adjacent to a charged surface in an electrolyte 
solution into motion. This classical subject of colloid sci- 
ence [l[ is currently experiencing a renaissance in micro- 
and nanofluidics 0, 13] , which raises fundamental ques- 
tion of how to pump and mix fluids at micron scales, 
where pressure-driven flows and incrtial instabilities are 
suppressed by viscosity. Electro-osmosis offers unique ad- 
vantages in this area of research and technologies, such as 
low hydrodynamic dispersion, no moving parts, electri- 
cal actuation and sensing, energy conversion and storage, 
and easy integration with microelectronics. 

Until recently, almost all studies of EO have assumed 
uniform surface charge and no-slip hydrodynamic bound- 
ary conditions at the surface. In such a situation the 
scalar electro-osmotic mobility Mi, which relates an ap- 
parent EO "slip" velocity Ui (outside of the thin double 
layer) to the tangential electric fleld Et is given by the 
classical Smoluchowski formula M 
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FIG. 1: (Color online) (a) Sketch of the superhydrophobic 
effective slippage effect on the EO flow. The real situation 
is approximated by a periodic cell of size L, with patterns 
of charges and flow boundary conditions (b) Illustration of 
tensorial EO response: 8 = ■k/2 corresponds to transverse, 
whereas 6 — to longitudinal stripes. 



For a charge density q2 of the slipping interface, simple 
arguments show that the EO mobility is given by 



X 



(1) 



where 77 is the viscosity of the solution, qi is the constant 
charge density of the no-slip surface, which can be related 
to the zeta potential across the diffuse (flowing) part of 
the double layer, Ci = qi/ne, where e is the permittiv- 
ity of the solution, and k = A^^^ is the inverse Debye 
screening length, that characterizes the thickness of the 
electrical Debye layer (EDL). 

Recent studies demonstrated the existence of a hydro- 
dynamic slip at hydrophobic smooth and homogeneous 
surface, which can be quantifled by the slip length h (the 
distance within the solid at which the flow profile extrap- 
olates to zero) 0] . The combination of the strategies 
of EO and hydrophobic slip, can yield enhanced EO fiow. 
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Since the EO flow ampliflcation scales as (1 -I- bn}, and 
b can be of the order of tens of nanometers [lO, U, 
12|, for typically nanometric Debye length an order of 
magnitude enhancement might be expected. 

It is now natural to assume that a massive amplifica- 
tion of EO fiow can be reached on super-hydrophobic sur- 
faces where effective, in general case tensorialj_slip length, 
boff, could be the order of several microns [ij, [ij, KLSf . 
The controlled generation of such flows is by no means 
obvious, since both the slip length and the electric charge 
distribution on a SH surface are inhomogeneous and of- 
ten anisotropic. Despite its fundamental and practical 
signiflcance EO flow over SH surfaces has received little 
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attention. Recently, [ig investigated EO flow past inho- 
mogcneously charged, flat SH surfaces in the case of thick 
channels {h L), thin EDL {kL ^1), and predicted 
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by using the Lorentz reciprocal theorem for the Stokes 
flow and by assuming perfect slip (6 = oo) at gas sec- 
tors. Here I is the unit tensor, and we keep notations, 
qi and 92, to characterize the surface charge density at 
the no-slip and slip regions, as above. This expression 
indicates negligible flow enhancement in case of an un- 
charged liquid-gas interface (which has been confirmed by 
later studies T3, 18 1 ) , and shows that surface anisotropy 
generally leads to a tensorial EO response. 

In this Letter, a general situation of EO flow past SH 
surfaces with patterns of arbitrary partial slip, is consid- 
ered (FiglD). Our focus is on the canonical EO geometry 
of a thick parallel-plate channel with a two-component 
(no-slip and slip) coarse texture, varying on scales larger 
than the EDL thickness. 

General theory.- To highlight the effect of 
anisotropy, we focus on an idealized, flat, periodic, 
charged, striped SH surface in the Cassie state, sketched 
in Fig. [U where the liquid-solid interface has no slip 
{bi = 0) and the liquid-gas interface has partial slip 
{b2 = b, < b < co). As a simple estimate, lubricat- 
ing gas sectors of height e with viscosity rjg much smaller 
than r] [l^ have a local slip length 62 ~ ^iv/Vg) ~ 50 e, 
which can reach tens of fj,m. Let then 0i and 02 ~ S/L 
be the area fractions of the solid and gas phases with 
ipi + 4'2 — 1- Pressure-driven flow past such stripes has 
been shown to depend on the direction of the flow, and 
the eigenvalues of the slip-length tensor [i^] read 21 1 
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These expressions depend strongly on a texture period L. 
When b/ L <Si 1 they predict the area- averaged isotropic 
slip length, 6^^" ~ 4>2b. When b/L '3> 1, expressions dU 
and (O take form 
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that coincides with results obtained for the perfect slip 
{b = 00) stripes [22 1. 



The EO mobility is represented by 2 x 2 matrices di- 
agonalizcd by a rotation [isj . By symmetry, the eigen- 
directions of M correspond to longitudinal (0 = 0) and 
transverse {9 = 7r/2) alignment with the applied electric 
field (Fig.[T]), so we need only to compute the eigenvalues, 
A/'l and M^, for these cases. 

We consider a semi-infinite electrolyte in the region 
y > Q above a patterned surface at y = subject to 
an electric field, Et = in the x direction. For 

nano-scale patterns (L < 1/xm), we can neglect convec- 
tion [Pe ^ 1 for a typical ionic diffusivity Z?), so that 
■>p{x,y,z) is independent of the fluid flow We also 

assume weak field {\Et\L <^ \tlj\) and weakly charged 
surface (|^| <^ ksT/lze) = 2b/ z mV) for a 2; : z elec- 
trolyte, so that -0 satisfies the Debye-Hiickel equation 
with a boundary condition of prescribed surface charge. 
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The fluid flow satisfies Stokes' equations with an electro- 
static body force 



0, 



(8) 



■qV u = Vp + SK il^EtX, V- 

with the boundary conditions at y = 

Ut = b{x, z)dyUt, ij ■ u = 0, (9) 

where Uf ~ ux + wz is the lateral, and v ^ y ■ u is 
normal to the surface velocities. We also neglect surface 
conduction (which tends to reduce EO flow) compared 



to bulk conduction {Du <^ 1) [23|. Far from the surface. 
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00, u approaches EO slip velocity U = — M ■ Et 
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For a longitudinal configuration only velocity component 
parallel to Et = EtX remains. In case of transverse to 
applied field stripes, normal velocity v ■ y does not vanish 
due to mass conservation condition in (|8]) , which can sig- 
nificantly modify the EO flow. Rigorous calculations [23[ 
allow one to find exact solutions for C/H'-'-, and thus ob- 
tain the eigenvalues of the EO mobility tensor: 
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where the effective slip lengths are given by Eqs.(j4]),([5]). 
The flow is thus anisotropic and there is a simple re- 
lationship between the EO mobility and hydrodynamic 
slip- length tensors 24 1 
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In the limit of 6/L ^ 1 the general expression transforms 
to Eq. ([3]). When &/X < 1 we get isotropic EO flow 



M = <t>iMi + (t>2M2 
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FIG. 2: Eigenvalues of normalized EO mobility. Solid curves 
represent longitudinal, and dashed - transverse alignment of 
stripes with electric field. Af "'^/Mi vs amplitude of the local 
slip b/L for (a) uncharged slip areas (§2 = 0,(1)2 = 0.5, HiL = 
10^), (b) uniform charge distribution (g2 = qi; 4>2 = 0.45, 
kL — 10'^) and (c) oppositely charged slip and no-slip areas 
(92 = -qi- 02 = 0.35, kL = 10^). M^^'-^/Mi plotted against 
(d) the fraction of gas sectors {b/L — 0.1, q2 — —qi, nL — 
10^). 



Discussion.- To demonstrate examples of very rich 
and unusual fluid behavior at the SH surface it is in- 
structive to consider some limiting cases of Eq. (|12p with 
different values of qi and q2- The results (shown in Fig. [2]) 
are somewhat remarkable. We see, in particular, that if 
gas area is uncharged, the EO flow related then only to 
the charge qi on the solid-liquid interface is generally in- 
hibited as compared with a homogeneous, solid no-slip 
surface with uniform charge density (see Fig. [H^a)) 
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We remark and stress that in contrast to common expec- 
tations the situation described by Eq.([T4|) corresponds 



to M" < A/-'-, i.e. the maximal directional mobility is 
attained in a transverse, and minimal - in longitudinal 
direction [25[. When 6/L <g; 1 we simply get M ~ 4>iMi. 
In other words, the (isotropic) EO mobility shows no 
manifestation of the slip, being equal to the surface av- 
eraged velocity generated by no-slip regions. This result 
coincides with expected for hydrophilic slip sectors. In 
the limit of 6/L ^ 1 this inhibition becomes negligibly 
small, and we obtain the simple result of [l6, 18 1, where 
EO mobility becomes equal to Mi regardless of the ori- 
entation or area fraction of the slipping stripes. These 
results suggest that although the absence of the screen- 
ing cloud near the gas region tends to inhibit the effective 
EO slip, the hydrodynamic slip acts to suppress this in- 
hibition. 



The situation is very different if the slipping interface 
carries some net charge, which is not an unreasonable as- 
sumption [i^ . To gain some insight into the possible EO 
flow enhancement, we consider first the case of uniform 
surface charge qi =92, where Eq. gives 



(15) 



M = • [I + Kbcft] 



which might be seen as a natural tensorial analog of 
Eq. ([2]). Fig. [2{b) includes theoretical results calculated 
with Eq. (jisp for a geometry of stripes, and is intended 
to demonstrate that the flow is truly anisotropic and can 
exhibit a large enhancement from effective hydrodynamic 
slip, possibly by an order of magnitudes. We stress that 
such an enhancement is possible even at a relatively low 
gas fraction, i.e. when boff is relatively small (but the am- 
plification factor, (I -I- Kboff), might be huge). Also note 
that in this situation Af H > M^, i.e. the fastest/slowest 
direction can correspond only to longitudinal/transverse 
stripes. 

An interesting scenario is expected for oppositely 
charged solid and gas sectors. If qi = —q2, then Eq. p2|) 
transforms to 
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which for b/L I simply gives M = Afi[(/)i — (/)2(1 + Kfe)]. 
The calculation results for this situation are presented in 
Fig.IDJc), and suggest a very rich fluid behavior. We see, 
in particular, that inhomogeneous surface charge can in- 
duce EO flow along and opposite to the field, depending 
on the fraction of the gas area as shown in Fig. [H^d) . Al- 
ready a very small fraction of the gas sectors would be 
enough to reverse the effective EO flow. Another strik- 
ing result is that electro- neutral surface (</>! 91 -1-^2^2 = 0) 
can generate extremely large EO slip. With our numeri- 
cal example this corresponds to 02 = 0.5. In other words, 
super-hydrophobic surface of average positive charge or 
even zero charge can induce an EO flow (different for 
longitudinal and transverse to applied field stripes) in 
the direction of the applied field as if it is uniformly and 
negatively charged. These findings are similar to those 
of y, I27 1 that the electrokinetic mobility depends on the 



charge distribution on the object, and not solely on its to- 
tal charge. However, in our case the flow is dramatically 
amplified due to hydrodynamic effective slip. 

These results may guide the design of SH surfaces 
for transverse electrokinetic flows in microfluidic de- 
vices . As we have shown above, effective EO mobility 
of anisotropic striped surfaces is generally tensorial, due 
to secondary flow transverse to the direction of the ap- 
plied electric field. Anisotropy {\Uz/Ux\) is maximized 
in certain direction 0max (as it is seen in FiglSK) and re- 
quires that (72/91 and nb are as large as possible. In a 
thick SH channel a transverse "plug" EO flow seems to 
be very fruitful direction compared to transverse hydro- 
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FIG. 3; The ratio of EO velocity components, \Uz/Ux\, (at 
(j)2 = 0.5, kL — 10^) as a function of (a) angle 9 at h/L = 10^ 
and (b) local slip length at the optimal angle 6. From top to 
bottom 172 /<Zi = 2, 1, 0.5 
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FIG. 4: (Color online) Streamlines of the EO flow computed 
at (1)2 = 0.35 and kL = 10^ for 92 /qi = -0.43 and 6 = 7r/2. 
The origin of coordinates coincides with the center of the gas 
region. From left to right the local slip length is 6/L = 0.1, 
0.5 and 10^ 



dynamic phenomena, where flow is "twisted" only near 
the wall j28l |. 

Another mixing mechanism is related to the formation 
of patterns of steady convective rolls on the scale propor- 
tional to the texture period (Fig. S]). This can happen 
in the situation of oppositely charged and transverse to 
applied field stripes. Fig. |4] illustrates the effect of the 
local slip on the morphology of the steady rolls forma- 
tion. We see, in particular, that increase of h leads first 
to appearance of additional convective patterns near no- 
slip areas, and then to a transition to a flow morphology, 
where recirculation of a fluid is observed only at the no- 
slip regions. This, in turn, induces the flow reversal (as 
in Fig.Hd)) 

Concluding remarks.— We have described EO on 
inhomogeneously charged and slipping anisotropic sur- 
faces. Our analysis provided the necessary tools to de- 
scribe a significant modification of EO phenomena on SH 
surfaces: to quantify the inhibition and enhancement of 
flow, the transition from its anisotropy to isotropy, onsets 
of convective rolls formation and a relevant flow reversal, 
which can generate a huge EO slip even in the situation 
of a zero mean charge. Our results may find numerous 
appHcations in microfluidic lab-on-a-chip devices. 
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